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Abstract 



We characterize the coexistence of itinerant ferromagnetism and spin- 
triplet superconductivity within a single mechanism involving local (Hund's 
rule) exchange among d electrons. The ratio of transition temperatures and 
the spin anisotropy of the superconducting gap is estimated for ZrZn2. The 
A phase is stable in very low applied and molecular fields, whereas the Al 
phases persists in higher fields. A small residual magnetic moment is present 
below the Stoner threshold in the superconducting phase. 
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The coexistence of weak itinerant ferromagnetism (WFM) and superconductivity (SC) 
has been recently discovered in UGe2 (Ref. 1), ZrZn 2 (Ref. 2) and URhGe (Ref. 3). The su- 
perconducting phase was encountered in the ferromagnetic phase; both of these phases seem 
to disappear with increasing pressure. Therefore, the superconductivity must be influenced 
by the ferromagnetism, particularly since the ratio of the Curie temperature (T c ) relative to 
the superconducting transition temperature (T s ) can exceed an order of magnitude. 

It is hard to imagine that the superconducting pairing in that situation involves a spin 
singlet, since the molecular field, e.g., in ZrZn 2 due to the exchange interaction, which is of 
the order of i H m = 150 T (fiBH m = 17 meV), exceeds by far the thermodynamic critical 
field H C <1T. 

In this paper we consider both of these types of ordering within a single mechanism - 
the Hund's rule exchange and we draw some universal conclusions from a relatively simple 
and testable model containing two microscopic parameters (apart from the density of states 
and its derivatives at the Fermi level). We follow some of the approximations of the original 
Bardeen, Cooper, and Schrieffer (BCS) theory, although we employ the linearized disper- 
sion relation and spin-split structure of quasiparticle states to account for the Fermi-liquid 
structure of the weak- ferromagnetic such as ZrZn 2 . 

The spin — triplet pairing in weakly ferromagnetic systems has been considered before 
being mediated by the exchange of longitudinal spin fluctuations! or being triggered by 
the electron-phonon interaction^! Quite recently&i, the question of the coexistence of fer- 
romagnetism with spin — singlet superconductivity has been reexamined within the mean- 
field theory. All the foregoing work is based on one-band (Hubbard or extended Hubbard) 
model, so the superconductivity can arise either from exchange of a paramagnon for repul- 
sive interactional, or as a result of local attractive interaction (negative — U model) B In 
our approach the pairing, induced by local Hund's' rule exchange, appears in the correlated 
and orbitally degenerate systems and together with the short-range Coulomb interaction is 
regarded as the source of itinerant magnetism in 3d and 4d metallic system. 0Th is interor- 
bital interaction remains local if the degenerate bands are not strongly hybridized. This is 
exactly what happens for ZrZn 2 , where the main contribution to the high density of states 
p(ep) at the Fermi energy from the two d bands, which do not hybridizei'S (see also the 
discussion below.) 

We start from the effective model proposed by us recently^ and extend it to consider 
explicitly both ferromagnetism and superconductivity. It is represented by the Hamiltonian 

H = £ ( E ke - V) n ki« + + n c (1) 

The quasiparticle band energy Eu t = — fi of quasimomentum k is labelled by the 
orbital index i = 1 and 2. The local interorbital and intraatomic Hund's' rule coupling can 
be represented in two equivalent ways asi 

H ex = -2JJ2(Sa ■ S i2 + -n u n 2i ) = -2J^A+,Am, (2) 

i im 

where Su and na are the spin and the particle-number operators for £-th orbital on site 
i, and af m are real-space spin-triplet creation operators on site i (Af x = af^a^p etc.). 7i c 
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represents the direct Coulomb interaction (the Hubbard term), as well as the interorbital 
Coulomb term 



(3) 



The last term does not influence the phases considered in the first order so it is dropped 
out. So, (1) represents, in our view, a correlated systemS, in which the local interactions 
determine the quantum instabilities. 

We introduce the combined Hartree-Fock-BCS approximation and the four dimensional 
Nambu-type representation. This means that the interaction part is rewritten first (up to a 
constant) in the following manner 



'Hex + T~(-c — 2 J ( < Ar im > A im -\- < A\ m > Ai m ~~ | < Ai m > 



(4) 



where S z =< Sfi + ^ti > is the magnetic moment per atom, < A im > is one of the three 
possible components of the superconducting gap parameters, and / = U + 2 J is the effective 
Stoner parameter. Subsequently, the Hamiltonian can be rewritten as a 4 x 4 matrix with 
creation operators fjj = (/fe 1T /fe u > f-k%v -f-ifel)' so that 



n = Y,f k Af k + Y,Ek2 + N 

k kcr 

where A is 4 x 4 matrix of the form 



im 2 + e 



< A-im > 

2J 



(5) 
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-E k2 + IS' , 
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A-i 





-E k2 - IS* J 



This matrix can be diagonalized analytically. The interesting cases are only those with the 
spin dependend gaps (i) Af > Aj_ and Ao = (which we call the anisotropic A phase); and 
(ii) Af 7^ and Aj = Ao = (which is called the AI phase). In case (i), the four eigenvalues 
are 



in -i 1/2 

X kal,2 = ^ E kl ~ E k2) T [4(^1 + E k2 - aISZ f + A l 



(6) 



The sign (=1=) corresponds to the labels (1,2) of \ ka and reflects hole and electron excitations, 
respectively. The spectrum is separated with respect to the spin orientation a — ±1 of the 
Cooper pair. The spectrum is fully gapped if both Af and Af are nonzero. We obtain 
a combination of the spin splitting and superconducting gap (the gap at the Fermi energy 

21 S z when the spin is flipped and Af or Af 



7*^2 + A 2 + 



{IS Z ) 2 + A 2 



E F = fi) is 2 (yV^; Tu t T Vl iu ; -r^ 

within the spin subband. What is probably more important, in Al state half of the spectrum 
remains gapless A^m = —E k2 + IS Z , = E kl + IS Z . Thus, there should be a substantial 
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linear specific-heat term present also in the superconducting Al state and this dependence 
should be distinguished from the T n dependence (n > 2) due to the gap zeros (the latter 
would require the interband hybridization and hence the dependence!!-] J — > Juu')- The 
linear specific heat 7 comes from the spin minority electrons, and since for a weak itinerant 
magnet IS Z <C Ep, there will be drop in the relative value of 7 of the order Pi/{p\ + Pi), 
where p a = p a (Ep) is the density of states in spin-a subband at the Fermi energy. 

The results obtained so far are general in the sense they are independent of a particular 
electronic structure. In the following we assume that the bands are the same, i.e. = 
Ey, 2 = Ey, so that Ey, — p ~ v F k, where Vp is the Fermi velocity. 

The Bogolyubov quasiparticle operators can also be easily calculated; then they areEl 

(7) 





with the coherence factors 



4] \ = ±( l + Wka 



V k 



y/2 \ 1 " K/\ ka 



(8) 



and the eigenvalue = \fc al = —^k a2 = [{ v <rk) 2 + A^.] 1 / 2 . Note that the spin dependent 
Fermi velocity is caused by the circumstance that we have spin-split bands in the ferromag- 
netic phase. The Hamiltonian has the diagonal form 7i = J2k a ^fccr(a^ ct^ a + P_f Sa /3_f £o .) + 
E , which is the BCS form with spin-dependent quasiparticle and gap energies. The last 
factor contributes to important differences with the BCS theory. Namely, the gap parameter 
at temperature T = is determined from the equation 



where p a = 12e Fa /W (ef a is the Fermi energy for quasiparticles with spin a and W is 
the effective width of the band, related to the Fermi velocity via W = (2Air 2 hv p/ilo) 1 ^ 3 , 
where Qq is the elementary- cell volume). The integration boundary k m is determined by the 
condition that the paired particles with spin a are present only within the spin-split region 
of the band, i.e., by the constraint v a k ma = aIS z . In effect, we obtain an estimate 

A - = 4 ^ exp ' (10) 

Analogously, we can estimate the critical temperature by selecting an equivalent but slightly 
different representation of the quasiparticle energies: Ey, — p — aIS z = hvpk — aIS z , where 
vf is the Fermi velocity in the paramagnetic phase. Then the condition for T$ (for A = Af) 
reduces to 



tanh I 
^ 

flVF 



hv F k-2,JS z 
2k B T s 



hvpk — JS Z 



(11) 



In effect, we obtain the estimate 
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2.26IS*exp[-—). (12) 



In both expressions for A CT and for Tg, the exponent contains the exchange integral, whereas 
the preexponential factor is multiplied by the magnetic moment. This means that spin-triplet 
superconductivity disappears together with ferromagnetism. This result is more general as 
it relies on the well founded notion that Cooper per of spin orientation a exist within the 
corresponding spin subbands if only A = (the B phase is not stable). 

To make the estimates explicit we have to relate the results for a weak itinerant fer- 
romagnet. It is easy to rederive those results in the present situation with T c ^> Tg. 
Namely, the Curie temperature is given by the expression 



T 

- 1 c 

7T 



V6 '-'^ 



P 



[(ip-iyip] 1 ' 3 , (is) 



with p' and p" being, respectively, the first and the second derivative of the density of 
states p(E) taken at Ep. Additionally, the magnetic moment at T = is given by S z = 
(l/2)[(/p-l)/5] 1 /2 ; W ith 



B = I* P - 




P 



(14) 



Taking the parabolic density of states corresponding to the linearized dispersion relation, 
we can determine the ratio T s /T c in an explicit manner 



Similarly, the gap ratio is 



l^mexpf-J-V (15) 



|l^exp((2//3J)( P 7p 2 )S«), (16) 



which in the limit Ip ~ 1 gives the ratio exceeding three. The ratio will grow rapidly with S z 
and the prediction ln(A|/A|) ~ S z could be tested experimentally. The Al superconducting 
state is reached fast with the increasing moment. 

To interpret our findings in the coexistence regime we can say that the A phase (with 
Aj 7^ 0) cannot appear either near the quantum critical point - the Stoner threshold, at 
which Ip = 1, or in the ferromagnetically saturated state. In the former case, the relatively 
small field can polarize totally the system (since the magnetic susceptibility x — Xo/(l ~~ Ip) 
is almost divergent). Under these circumstances, the bound state with Aj 7^ cannot be 
formed, since the polarized surrounding will respond to the presence of the second electron 
strongly. On the other hand, if the magnetic moment in the ferromagnetic state is almost 
saturated, then the bound state with Aj_ cannot be formed as the system is rigid. Only 
in between those two limiting situations, but rather on the weak-ferromagnetic side, the 
coherent A state can be realized. Otherwise, the Al state will become stable. 
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To gain a quantitative insight into the nature of A and Al states and their coexistence 
with the spin polarized state, we have also performed the analysis in the paramagnetic 
state, i.e., for pi < 1, but in the applied magnetic field B ^ 0. To amplify numerically the 
effects discussed, we have put I = J and choose a constant density of states . The phase 
diagram as a function of magnetic field is displayed in Fig. 1. Note again that the phase A 
disappears in the vicinity of the Stoner threshold marked by the vertical dashed line. On 
the contrary, the phase Al continues towards the Stoner boundary. In Fig. 2, we provide 
the field dependence of the ground state energy away from the Stoner threshold for normal, 
A, and Al phases. The Al phase is the most stable in a high applied field and can coexist 
with a saturated ferromagnetism. The same happens for the superfluid 3 He. This means 
that the superconducting coherence length for Al phase becomes unbound when the system 
reaches the saturated state. 

The interesting feature of our results is that the presence of the superconductivity indices 
a magnetic moment even below the Stoner threshold. The magnitude of the moment below 
the Stoner threshold ( J/W = 0.25 in this case) is displayed in Fig. 3. Obviously, the effect is 
enhanced by the choice of the density of states selected and particularly by the assumption 
I = J , but it is worth mentioning. It suggests that in the presence of the superconducting 
state makes the Stoner critical point a hidden one. This point will be elaborated further 
elsewhere. 

In summary, we have considered the most natural model for the coexistence of ferro- 
magnetism and spin-triplet superconductivity, both induced by a single mechanism - the 
local ferromagnetic exchange in the orbitally degenerate systems. While the anisotropic 
A-state seems to be stable in weakly polarized and paramagnetic systems at low fields, the 
Al state may coexist even closer to the quantum critical point. The paired state induces 
a small magnetic moment even below the Stoner threshold. It should be interesting to ex- 
tend these results to the hybridized systems^, when fc-dependent (p- or d-wave) spin-triplet 
superconductivity will appear. 

This work was supported by the KBN Grant No. 2P03B 092 18, as well as by NSF Grant 
No. DMR 96-12130. 
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Figure Captions 



Fig.l. Phase diagram specifying the stable superconducting states A and Al (see main 
text). The dotted line specifies the onset of the magnetically saturated state. The vertical 
dot-dashed line marks the Stoner threshold. 



Fig. 2. Ground-state versus applied magnetic field B. At B = the equal-spin pairing 
state (A-phase) is stable, whereas in high field Al phase (with the spins parallel) prevails. 
The system is below the Stoner threshold. 

Fig. 3. Magnetic moment per orbital (upper panel) and the superconducting gaps (lower 
panel) versus applied field (below the Stoner threshold). The inset shows the chemical 
potential vs. B. The horizontal arrow indicates a small residual magnetic moment in the 
B —> limit below the Stoner threshold. 
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FIG. 1. 
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FIG. 2. 
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